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A potential explanation [U. D. Jentschura, Phys. Rev. A 88, 062514 (2013)] of the proton radius 
puzzle originating from the non-perturbative lepton-pair content of the proton is studied. Well- 
defined quantities that depend on this lepton-pair content are evaluated. Each is found to be of 
the order of 10(^) 2 , so that we find such a lepton-pair content exists in the proton. However, we 
argue that this relatively large result and general features of loop diagrams rule out the possibility 
of lepton-pair content as an explanation of the proton radius puzzle. The contributions of a class 
of potential explanations of the proton radius puzzle (for which the dependence on the pp relative 
distance is as contact interaction) are shown to be increase very rapidly with atomic number. 

I. INTRODUCTION 

Recent high precision experimental studies of muonic hydrogen urn obtain a value of the proton radius that is 
about 4% smaller than that obtained from ordinary electronic hydrogen. The problem of understanding this difference 
has become known as the proton radius puzzle, and has generated a vast array of possible solutions, see the review 
[5j. One of the novel suggested solutions 0 is that a non-perturbative feature of the proton’s structure, namely the 
possible presence of light sea fermions as constituent components of the proton, could account for the difference in 
the extracted radii. In particular, it is argued that the assumption that the presence of 2.1 x 10 -7 light sea positrons 
per quark, leads to an an extra term in the electron-proton versus muon-proton interaction, which has the right sign 
and magnitude to explain the proton radius puzzle . 

The basic idea is that a bound electron may annihilate with a positron that is part of the non-perturbative e - e + 
cloud of the proton. The annihilation leads to a virtual photon, which in turn decays to a bound electron and a 
positron that is also part of the e _ e + cloud of the proton. See Fig. [l] The term non-perturbative here refers to a 
component of the proton Fock-space wave function that can be seen at small values of momentum transfer. If one 
could take a snapshot of a proton in isolation one would see electron-positron or muon-anti-muon pairs pop in and 
out of existence. This effect is therefore different from the generation of pairs by evolution in momentum transfer 
that is akin to the source of the qq sea of perturbative QCD. 



FIG. 1: (color online) Typical Feynman diagram of Ref. j4] illustrating the virtual annihilation of a bound electron with a 
“light sea lepton” (positron) inside the proton. The up ( u ) and down (d) quarks, which carry non-integer charge numbers, 
interact electromagnetically. The vertical dashed lines indicate the lepton-sea component of the proton wave function. The 
light sea lepton that annihilates with the bound electron. An example of a non-perturbative quantum chromodynamic (QCD) 
interaction via a blue-antigreen gluon also is indicated in the figure. 

A natural question to ask is whether or not this diagram is part of a contribution that is already included. In 
particular, a look at Fig [I] might lead one to conclude that the intermediate baryonic state is that of an excited 
nucleon. If so, this diagram would be a particular time-ordering of the proton polarizability contribution to the two- 
photon exchange diagram. However, if the lepton pair is a specific Fock-space component of the complete proton wave 
function (including QED effects) one may argue j3] that the intermediate state is part of the proton wave function, 
and therefore not part of the proton polarizability contribution. Here we accept this argument and seek to determine 
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its logical consequences. 

This acceptance comes with severe difficulties. If one adopts this Fock-space approach, one must recognize that 
computing the contribution of a given component may only correspond to computing a particular time ordering of 
a Feynman diagram, with the consequent obligation to find any remaining time ordering terms of the same order in 
a. To be specific, the term of Fig. 1 corresponds to a particular time ordering of the graphs that give the radiative 
corrections to the general two photon exchange diagrams. See for example, Fig.pJ In this figure, the intermediate blob 
is meant to contain the proton and all of its excited states. Thus this graph contains the time-ordering that is included 
in Fig. 1. One may immediately estimate the size of the contribution of this term to the energy to be of order a/n 
times the polarizability correction, and is therefore negligible. The explicit calculations for muonic hydrogen are in 
accord with this estimate [3 El- For electronic hydrogen the contribution is again of order a/ir times the polarizability 
correction, which is relatively much less important for electronic hydrogen than for muonic hydrogenjT]. This means 
that a correct evaluation of the contribution of the effects of the proton’s lepton-pair content must give a negligible 
result. 


e 


P 


FIG. 2: (color online) Radiative correction to two photon exchange. Including the full gauge invariant set, including crossing 
the photons arising from the proton (p) and all of the one-loop self-energy terms on the electron(e), is implied. 

We pursue the idea of Ref J4] even though the result of any correct calculation will be that the effect is negligible. 
This is because investigating the lepton-pair content of the nucleon is interesting in its own right and because it useful 
to see how the correct result arises from the Fock-space idea. The argument of Ref [3] starts by considering the photon 
annihilation term occurring in positronium. This term leads to the effective interaction [§]: 

6 H =^(3 + <T + -cT_)6(r). (1) 

This Hamiltonian gives a nonzero interaction of the bound electron and the light sea positron if their spins add up to 
one. Ref. [3] assumes that the electron-positron pairs within the proton are not polarized and replaces <7+ • c?_ —> 0 
after averaging over the polarizations of the sea leptons. Further it is argued [4j that for atomic (electronic) hydrogen, 
the additional interaction of the electron with the proton due to the annihilation channel takes the form 

H ann = e p ^*6(r), (2) 

where e p measures the amount of electron-positron pairs within the proton. For muonic hydrogen, the effect is 
expected to vanish because the dominant contribution to the sea leptons comes from the lightest leptons, namely, 
electron-positron pairs and thus the annihilation channel is not available. Ref. [1] finds that a value of e p = 2.1 x 10 ' 
is sufficient to account for the different values of the extracted radii. Note that the 1 /m 2 e dependence occurs as the 
result of assuming that both the electron and positron are at rest. 

The hypothesis of sea leptons and the use of Eq. (pi) raises a number of interesting questions. Is the quantity e p 
well defined? If so, what is its likely range of values? Any positron in the non-perturbative nucleon sea is not likely to 
be at rest, so one could also ask if the Eq. (p]) is applicable at all. Does the term of Eq. p]) really exist? Our purpose 
here is to investigate these questions. 

Here is an outline of the remainder of this paper. The difficulties in computing and defining e p are discussed 
in Sect. II. Observables that depend on lepton-pair content and can be computed in a gauge invariant manner are 
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discussed in Sects.Ill and IV. The results of these two sections are analyzed and used to understand more detailed 
loop calculations Sec. V. The nuclear dependence of all models in which the contribution to the Lamb shift arises 
from an interaction containing a Dirac delta function in the lepton-nucleon separation is studied in Sec. VI. Some 
concluding remarks are made in Sect. VII. 


II. ATTEMPTING TO COUNT POSITRONS IN THE PROTON 

We might proceed to count positrons by defining an anti-lepton number current density operator, JF, 

= $n ( 3 ) 

and take its expectation value in the physical proton wave function. The operator Jk is defined to act only on 
anti-leptons (/), so its O’th component is the anti-lepton density. We may take the expectation value of this operator 
in the proton by assuming that the lepton pair arises from the interactions on a single quark, as in Fig. 2. In that 
diagram j represents the incoming momentum, p and p' = p + j are respectively the incoming and outgoing quark 
momenta, m and m q are respectively the lepton mass and the constituent quark mass, q is the momentum flowing 
into the inner loop, and k is the momentum on one branch of the inner loop. 

It is necessary to discuss why we focus on a single quark. In general the momentum in the loop q is very large. Any 
operator that removes a large momentum from one quark and gives a large momentum to another quark, suffers a 
vastly reduces matrix element in the proton wave function because each quark can only support a momentum of the 
order of the inverse radius of the proton, or about 200 MeV/c. There are significant implications to this single-quark 
dominance see below in Sect. eh 



(a) (b) 


FIG. 3: (color online) The solid line represents a single quark, (a) Diagram for counting anti-leptons. The notation X with the 
arrow coming into it represents the operator bringing in a four-momentum j. (b) Lepton-pair contribution to the vertex 
function. 

One immediate problem with the amplitude of Fig. [3] is that the integrations over the virtual lepton momenta 
contain an ultraviolet divergence. This means that an unknown counter term is needed to determine a value, and 
predictive power is lost. Furthermore, the use of of Eq. |3j) does not respect current conservation, because Jj does 
not act on all of the charges. Thus any result would be gauge-dependent. These two drawbacks lead us to conclude 
that the short answer to the question, “is the quantity e p well defined?”, is simply no. Thus the validity of Eq. (|2j) , 
which contains a factor of e p , is questionable. 

However, one may compute matrix elements that depend on the pair content that are free of ultraviolet divergences 
and are gauge invariant. This is the procedure adopted here. 


III. ANOMALOUS MAGNETIC MOMENT 

The goal of this section is to find a simple gauge invariant quantity free of ultraviolet and infrared divergent terms, 
that depends on the lepton-pair content of the proton and to evaluate that quantity. Such an example of a leading- 
order effect of pairs is shown in Fig. 3b, in which the X represents the usual electromagnetic current operator. The 
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intermediate electron line exists only if the lepton-pair exists, and computing only the contribution to the magnetic 
moment insures that gauge invariance is maintained without encountering an ultraviolet or infrared divergence. The 
effects of this diagram have been evaluated long ago, so our purpose is merely to find an illustrative example. 

The effect of pairs is accounted for by dressing the photon propagator using [B- 


Wfiiy 


q 


2 


-ig^ i 
q 2 l~Il 2 (q 2 ) 



(i + n 2 (<z 2 )), 


(4) 


for a virtual photon of four-momemtunr q , with 


n 2 (g 2 ) = 


-2 a 


dzz(l — z) log 


mf 


mf — z( 1 — z)q 2 ’ 


where mi is the lepton mass and the subscript 2 stands for second order in e. It is useful to rewrite this term as 


n 2 (g 2 ) = 


-2 a 


r \ n , 2 f 1 d\ 
dz z(l — z)q~ — 


1 


A mf — q 2 ’ 


where 


mf = 


mf 


z(1 — z )A 


(5) 


( 6 ) 


(7) 


The use of Eq. ([ 6 ]) allows the evaluation of a Feynman diagram involving II 2 via the usual technique of combining 
denominators. 

The order e 4 correction term to the vertex function, 6T^(q) is given by using Eq. Q in the standard expression |Bj 
for the e 2 term: 


5T fl (q)=2ie 2 


^ 4 ^ u{jp') (jkj 1 */;' + m 2 7 ^ — 2 m(k + k') /i ') u(p ) 
(27t) 4 ((p — fc ) 2 + ie)(k' 2 — m 2 + ie)(k 2 — m 2 + ie)) 


n 2 ((p-fc ) 2 


( 8 ) 


where p is the initial quark momentum, p' = p + q, k! = k + q and m is the quark mass, which we take at a constituent 
value of one-third of the mass of a proton. The use of Eq. © in Eq. Q gives the result: 


, 2a 


dX 


<5T M (g) = 2*e 2 — / dz z(l — z) / — 


u{p') + to 2 7 m — 2 m(k + k ')^ u(p) 


0 A J (27t) 4 (Jj) — k) 2 — mf + ie){k' 2 — m 2 + ie)(k 2 — m 2 + ie )) 


(9) 


The above correction term contributes to both the Dirac and Pauli form factors of the quark. Evaluating the Dirac 
form factor requires treatments of infrared and ultraviolet divergences that are not present in the Pauli form factor. 
We therefore compute only the contribution to the Pauli form factor. Then we proceed by combining denominators, 
and integrating over the four-momentum and two of the Feynman parameters and A. 

The relevant momentum transfer for atomic physics is the inverse of the Bohr radius, this is much, much smaller 
than the inverse of the quark mass. Therefore we need only evaluate 5T^(q) at q 2 = 0. We find the following result 


ST^q 2 = 0 ) = 2 ^) 2 u(p')^u(p)I{m 2 /m 2 ) 


I{m 2 /m 2 ) = dz z(l - z) f, 1 dz' log (l + 

Numerical evaluation leads to the result 

/(electron) = 1.7, /(muon) = 0.17, 


( 10 ) 

( 11 ) 

( 12 ) 


and the virtual e + e _ pair contribution to the magnetic moment of the quark is given by 3.4 (^) 2 or 1.8 x 10~ 5 . This 
number is about 100 times larger than the value of e p used to account for the proton radius puzzle. Moreover, if the 
value of the lepton mass in the loop were to dominate the value we would expect a muon to electron ratio of about 
1/200 instead of the ratio ~ 1/10 of the values of Eq. (20 1 . 

It is useful to attempt to relate the results of Eq. (20 1 to access the value of e p , even though this can only be done 
heuristically. In non-relativistic quantum mechanics the contribution of a component n to the expectation value of an 
operator O is written as (O) = f d 3 rifi*(r)Oi/> n (r). This can be interpreted as f d 3 r\4> n (r)\ 2 x <D = P n O. where O is 
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the average of O in the component n. If this operator is written in natural dimensionless units, such as a coefficient 
multiplying the factor u(p') zc7 2m giy u(p), then O can be expected to be of the order of unity. In that case the numerical 
coefficient, here 1.8 x 10~ 5 can be expected to be of the order of the probability for the electron-postiron pair to exist. 
Thus, we state, very roughly, that 

e p ~ 1.8 x 1(T 5 . (13) 

Note that the dimensionless factor, /, multiplying zt(p') J" u(p) depends on the ratio of the quark to lepton 
masses. This is not surprising-at q 2 = 0 the only parameters with dimension are m, mi and the only way to make a 
dimensionless number is a dependence on the ratio. The consequences of this are discussed below in Sect. |V| 


IV. AXIAL COUPLING 


We seek another example of a divergent-free matrix element that depends on the lepton-pair content of the proton. 
Consider the axial coupling 7 ^ 7 ® operator as an insertion in Fig.2b. The resulting term is defined as 5 A with 


5 A» = ie 2 


d 4 k Yp'Yp (V + m)YY(t + m)) YYp) 
( 2 tt ) 4 ((p — k) 2 + ie)(k ' 2 — m 2 + ie)(k 2 — m 2 + ie)) 


n 2 ((p-fc) 2 ), 


Using parity conservation and time reversal invariance tells us that 8A 44 takes the form 


SA 41 = u(p') 


G A (q 2 ) YY + 7 5 + 


i(p). 


(14) 


(15) 


The term G Al as computed from Eq. (14) will contain an ultraviolet divergence, and are determined after a renor¬ 


malization procedure. The terms GV, G A are free of such problems, so we will examine only those terms. We again 
take q 2 = 0 . 

We proceed by using Eq. © in Eq. ([l4|) to find 


, 2a 


dX 


SA tJ, =ie "— / dzz(l — z) / — 


d A k u( ,P'hp (y + m)Yl b { t + m)j 7 P u(p) 


o A J (27t) 4 ((p — k) 2 — hi 2 + ie){k' 2 — m 2 + ie)(k 2 — m 2 + ie)) 


(16) 


We proceed by combining the propagators, shifting the origin of integration, integrating over the shifted four- 
momentum, integrating over one of the Feynman parameters, taking q 2 = 0 and keeping only the contribution to 
Gr,Gp. The result is 


1 fa -\2 ^ f 1 ^ r 1 dX r 1 r 1_z ' i 2 + z ' i- 2 V z '+Z (l-z')) 


8A^ P (q 2 = 0) = h_(^) 2 mf 0 dzz(l-z)f 0 f f 0 dz'f n ~ dy- 


(l—z') 2 m 2 -\-z'rhf+ie 


(17) 


The integration over y gives a factor of 1 — z' to to term proportional to q l>J Y ■ but leads to a cancellation in the term 
proportional to ia v ‘ v q v Y■ Thus the Gp term vanishes at q 2 = 0. This cancellation does not occur for other values of 
q 2 . The integration over A is performed to give the result: 


,L4£( 9 2 =0) = Hf) 2 q -^Am 2 /m 2 ) 

J(m 2 /m 2 ) = dzz( 1 z) dz' « 2+ £ 1 +* , » log[l + ] 


w Jo 

The remaining integrals are handled numerically. The result is that 

J(q 2 = 0, electron) = 18.3, J(q 2 = 0,muon) = 0.593, 


(18) 

(19) 

( 20 ) 


We again see that the effect of e + e pairs is much larger than the 2 xlO ' that is supposed to enter in the proton 
radius puzzle. Using the logic of the previous section leads to the approximate relation 


5 x 10 


-5 


( 21 ) 


This has about 250 times larger than the value used in Ref. [1]. The numbers displayed in Eq. (13) and Eq. (21) are 
of the same order of magnitude, which is all that can be expected from the quantitative approach used here. 
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V. ASSESSMENT 


We have provided arguments that the lepton pair content is an order of magnitude larger than the value of e p 
needed to account for the proton radius puzzle. One might expect that this strengthens the case for the explanation 
of Ref. [3] . 

That this is not so can be understood by examining the validity of Eq. ([2]) , derived by assuming that the electron 
and positron annihilate when both are at rest. Firstly, if this equation is correct, and the lepton-pair probabilities 
are as large as obtained in this paper, then the computed effect would be between 100 and 250 times times too big., 
causing a disagreement with experiment that would rule out the explanation of Ref. [1] . But it is necessary to assess 
the lepton mass dependence that arises from evaluating the loop diagrams. This is in both examples approximately a 
logarithm of the ratio m 2 /mf, as shown in Fig. [4j This means that a dependence as the inverse square of the lepton 
mass will not result from evaluating any loop diagram. 


l(R)/Log(R) 



/ 

/ 

/ 


5.R=m 2 /m? 


FIG. 4: (color online) Dependence on the ratio ra 2 /ra 2 . The functions /(ra 2 /ra 2 )/log(ra 2 /raf ),Eq. (11) (blue, solid) and 
J(m 2 /m 2 )/ log(m 2 /m 2 ), Eq. (191 (red, dashed) are displayed. 


Moreover, we may use dimensional analysis to argue that a complete evaluation of the effect depicted in Fig. 1 
would lead to a 1/m 2 dependence. This is because the inverse factor of appearing in Eq. Q arises from the 
propagator of the virtual photon that is produced by the annihilation. In the diagram this factor would be (p~ +p+) 2 . 
The momentum p~ of the external electron is given almost entirely by the electron mass, but the four-momentum of 
the space-like virtual positron, p\ i", is dominated by its three-momentum. This in turn is governed in loop integrals 
by the largest mass, which is that of the quark, m. In particular, a correct and complete evaluation would lead one 
to obtain 

7 TCL 

SH ann oc (3 + er + • er_) <5(r) (22) 

instead of Eq. 0. Since m ss 600 rn e the effect is about four hundred thousand times smaller than what is needed 
and therefore can be said to vanish. 

We have used the term “a full and complete evaluation”. What would be needed for that? The minimum requirement 
is gauge invariance. That the diagram of Fig. 1 can not satisfy gauge invariance by itself is apparent because there 
are many other diagrams of the same order. In particular, there is a diagram corresponding to crossing of the photons 
emitted by the u and d quarks of Fig 1. Terms with both crossed and uncrossed photons are needed to satisfy gauge 
invariance in Compton scattering. This means that the lepton pair content can not be evaluated without also including 
the effects of the 2ye _ e + component. One simply needs to compute the complete gauge invariant set of diagrams 
corresponding to those implied in Fig 2. This has been done for the case of muonic hydrogen. Given the insensitivity 
of the results to the value of the lepton mass, we expect that for electronic hydrogen the lepton-pair effect will be of 
order of a/n of the proton polarizability correction and therefore negligible. 


VI. NUCLEAR DEPENDENCE 

Models such as the contact interactions of Eq. 0 which behave as a delta function in the separation between the 
lepton and nucleon contain very specific predictions for the nuclear dependence of the Lamb shift. 
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The first example we consider is the model of Ref. [2|. To make a prediction for the nucleus, one needs to know 
the contribution of the neutrons. Using the single-quark dominance idea of Sect II gives a specific result obtained 
by considering the factors of the square of the quark charge that would appear in the calcuation. A proton contains 
two up quarks and one down quark, with a resulting quark charge squared factor of 2(2/3) 1 2 e 2 + l(l/3) 2 e 2 = e 2 . For 
a neutron one would have l(2/3) 2 e 2 + 2(l/3) 2 e 2 = 2/3e 2 . Thus the neutron contribution would be 2/3 that of the 
proton. As result, if one considers the Lamb shift in the electron-deuteron atom, the effect would 5/3 as large as for 
a proton. Such an effect would contradict the existing good agreement between theory and experiment [Ill- 

More generally, suppose the contribution of a proton to the Lamb shift is E p (0.3meV) to resolve the proton radius 
puzzle) and that of the neutron is E n . Then for a nucleus with A nucleons and Z protons, we find 

/ i + St \ 3 / „ 2 x / i _(_ I'hi \ 3 

Ea = —Z 3 (ZE p + NE n ) 1 - -4 ) ~ Z 3 (ZE p + NE n ), (23) 

\ ' Am p ) V a ii J yd” Am p ) 

where a p is the muon Bohr radius (>100 times larger than nuclear radius, Ra )■ The meaning of Eq. ( |23[ ) is that the 
contributions of such contact interactions increase very rapidly with atomic number. 

In particular, the prediction of Ref. 0](with E n = 2/3 E p ) for 4 He is a Lamb shift that is (1-27) 8 (2)5/3« 10 meV, 
a huge number. The expression Eq. ( |23[ ) applies to all models in which the contribution to the Lamb shift enters 
as a delta function (or of very short range) in the lepton-nucleon coordinate, including (9} HO]. In Ref. [9] , which 
concerns polarizability corrections, the neutron contribution, E n could vanish, so that the contribution for 4 He would 
be 20 E p =6 meV. In Ref. [TO], which concerns a gravitational effect, E n = E p . so the prediction for 4 He would be 
40E p =12 meV. These various predictions will be tested in an upcoming experiment m- It is also worth mentioning 
the MUon proton Scattering Experiment (MUSE) p3j, a simultaneous measurement of p + p and e + p scattering and 
also a simultaneous measurement of p~ p and e~ p scattering that is particularly sensitive to the presence of contact 
interactions [9] . 


VII. SUMMARY 

The work presented here supports the idea of the existence of a non-perturbative lepton-pair content of the proton. 
Such components are not forbidden by any symmetry principle and therefore must appear. However, the presented 
calculations show that such a content is not a candidate to explain the proton radius puzzle. This is because 
computation of the necessary loop effects is cannot yield an effective Hamiltonian of the strength and form of Eq. (|2j) . 
The 1/m/ behavior of that equation in Ref. [2] is necessary to obtain the needed magnitude of the separate electron 
and muon Lamb shifts. Instead, loop calculations are expected to lead to a a dependence of 1/m 2 , with m the 
constituent quark mass. This means that the effect of Ref. [J is expected to be entirely negligible. This is in accord 
with the estimate: (cx/n) times the polarizability correction that is obtained from evaluating the relevant Feynman 
diagrams. 

More generally: it can be said that the proton does have Fock -space components containing lepton pairs. However, 
the simplest and most reliable method of treating such pairs is to compute gauge-invariant sets of Feynman diagrams 
using QED perturbation theory. 
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